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^ [ Abstract 

o , 

O ■ The Henon-Heiles system in the general form is studied. In a nonintegrable 

^ i case new solutions have been found as formal Laurent series, depending on 

^1 three parameters. One of parameters determines a location of the singularity 

^ ■ point, other parameters determine coefficients of the Laurent series. For some 

^ '. values of these two parameters the obtained Laurent series coincide with the 

^ ' Laurent series of the known exact solutions. 

m ; 
(N ■ 

> 1 The Henon— Heiles Hamiltonian 

m ■ 

o : 

Q , Let us consider a three-dimensional galaxy with an axial-symmetric and time-independent 

2 I potential function. The equations of galactic motion admit two well-known integrals: 

^ ' energy and angular momentum. If we know also the third integral of motion, then we can 

O ■ solve the motion equations by the method of quadratures. Due to the symmetry of the 

^ . potential the considered system is equivalent to two-dimensional one. However, for many 

^! polynomial potentials the obtained system has not the second integral as a polynomial 

^ [ function. 

^ ' In the 1960s numerical |jl| and asymptotic methods 0, have been developed to show 

H ■ either existence or absence of the third integral for some polynomial potentials. To answer 

. the question about the existence of the third integral Henon and Heiles [|l| considered 

^ ! the behavior of numerically integrated trajectories. They wrote []I|: "In order to have 

^ \ more freedom of experimentation, we forgot momentarily the astronomical origin of the 

■ ■ ■ ' problem and consider its general form: does an axisymmetrical potential admit a third 
isolating integral of motion ?" . They have proposed the following Hamiltonian: 

H = l (x? + + x^ + y^) + x^y - ^y\ (1) 

because: (a) it is analytically simple; this makes the numerical computations of tra- 
jectories easy; (b) at the same time, it is sufficiently complicated to give trajectories, 
which are far from trivial. Indeed, Henon and Heiles found that for low energies this sys- 
tem appeared to be integrable, in so much as trajectories (numerically integrated) always 
lay on well-defined two-dimensional surfaces. On the other hand, they also obtained that 
for high energies many of these integral surfaces were destroyed and that phase space 
acquired large ergodic regions. The Henon-Heiles system became a paradigm of chaotic 
Hamiltonian dynamics. 
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H 



i (x^ + + Xx^ + y^) + x'y - (1) 




Let us consider the Henon— Heiles Hamiltonian in the general form: 

C 

where C and A are numerical parameters. 

Investigations of the generalized Henon-Heiles system: 

—Ax — 2xy, 

(2) 

- + Cy"^, 

can be separated on the following ways: 

1. Numerical analysis |^ showed, that in the original case (A = 1, C = 1) 
singular points of solutions of the motion equations group in self-similar spirals. It turns 
out extremely difficult distributions of singularities, forming a boundary, across which the 
solutions can not be analytically continued. Numerical investigations of the generalized 
Henon-Heiles system are continued up to now P, 0. 

2. The procedure for transformation the Hamiltonian to a normal form and for con- 
struction the second independent integral in the form of formal power series in the phase 
variables x, Xt, y yt (Gustavson integral) has been realized for the Henon-Heiles system 
both in the original (A = 1, C = 1) (see also 0) and in the general forms 0. Using the 
Bruno algorithm [|TU], |TT| V.F. Edneral has constructed the Poincare-Dulac normal form 
and found (provided that all phase variables are small) local families of periodic solutions 
of the Henon-Heiles system both in the original [Q and in the general [0 forms. 

3. The singularities at the fixed points in phase space are locally analyzed via 
normal form theory, whereas the singularities in the complex (time) plane are studied 
by the Painleve analysis. Using this analysis three integrable cases of the generalized 
Henon-Heiles system have been found: 



(i) 


C = 


-1, 


A = 


1, 




C = 


-6, 


A is 


an 


in) 


c = 


-16, 


A = 


1 

16' 



In contradiction to the case (i) the cases (ii) and (iii) are nontrivial, so the integrability 
of these cases had to be proved additionally. In the 1980's the required second integrals 
were constructed [14-18]. For integrable cases of the Henon-Heiles system the Backlund 
transformations |l^ and the Lax representations [^, ^] have been found. The three 
integrable cases of the Henon-Heiles system correspond precisely to the stationary fiows 
of the only three integrable cases of fifth-order polynomial nonlinear evolution equations 
of scale weight 7 (respectively the Sawada-Kotega, the fifth-order Korteweg-de Vries and 



the Kaup-Kupershmidt equations) |2T|, |22|. The general solutions of the Henon-Heiles 
system in integrable cases are known . 

4. The Henon-Heiles system as a system of two second order ODEs is equivalent to 
the fourth order equation|: 

20(7 

yutt = (2C - ^)yuy - (4A + l)yu + 2(C + l)yl + ^y^ + (4CA - Q)y^ - Xy - AH. (3) 



^For given y(t) the function is a solution of a linear equation. System (2) is invariant to exchange 
X to —X. 
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To find a special solution of this equation one can assume that y satisfies some more 
simple equation. For example, it is well known that the Henon-Heiles system and, hence, 
equation (3) have solutions in terms of the Weierstrass elliptic functions satisfying the 
first-order differential equation: 

yl = Ay^ + By^ + Cy + V, (4) 

where A, B, C and V are some constants. 

E.l. Timoshkova P3 generalized equation (4): 



y^ = Ay' + By' + Cy + V + Gy'/' + Sy'/' (5) 

and found new one-parameter sets of solutions of the Henon-Heiles system in noninte- 
grable cases (C = — | or C = — y, A is an arbitrary number). 

In the present paper I use the Painleve method to find asymptotic solutions of the 
Henon-Heiles system at C = — y. 



2 The Painleve property 

Let us formulate the Painleve property for ODE's. Solutions of a system of ODE's are 
regarded as analytic functions, may be with isolated singular points ||2^, A singular 
point of a solution is said critical (as opposed to noncritical) if the solution is multi- 
valued (single-valued) in its neighborhood and movable if its location depends on initial 
conditions^]. 

Definition. A system of ODE's has the Painleve property if its general 



solution has no movable critical singular point ^7 . 



An arbitrary solution of such system is single- valued in the neighborhood of its singular 
point to and can be expressed as a Laurent series with a finite number of terms with 
negative powers of t — to- 

A Hamiltonian system in a 2s-dimensional phase space is called completely integrable 
if it possesses s independent integrals which commute with respect to the associated 
Poisson bracket. When this is the case, the equations of motion are (in principal, at 
least) separable and solutions can be obtained by the method of quadratures. Since 
the work of S.V. Kovalevskaya ||2^ (see also [^) on the motion of a heavy rigid body 



about a fixed point, the Painleve property has been proposed as a criterion for complete 
integrability [30-32] . If the system misses the Painleve property (has complex or irrational 
"resonances"), then the system cannot be "algebraically integrable" (see also and 



references there in). N. Ercolani and E.D. Siggia ||35| , pq| advance arguments as to why 
the Painleve test works, i.e. they showed how to exploit the singular analysis to yield 
the integrals. They proved a theorem which demonstrates that the singularity analysis 
provides bounds on the degrees of polynomial integrals for a large class of separable 
systems. 



^Solutions of a system with a time-independed Hamiltonian can have only movable singularities. 
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3 The Painleve test 



3.1 Various algorithms of the Painleve test 

The Painleve test is any algorithm designed to determine necessary conditions for a dif- 
ferential equation to have the Painleve property. The original algorithm developed by 
Painleve |27| is known as the a-method. The method of S.V. Kovalevskaya is not as 



general as the a-method but is more simple than it is. 

In 1980, motivated by the work of S.V. Kovalevskaya p8[, M.J. Ablowitz, A. Ramani 



and H. Segur developed a new algorithm of the Painleve test for ODE's 0. Using this 



algorithm one can determine whether a system of ODE's admits movable branch points, 
either algebraic or logarithmic. This algorithm can be used not only to isolate values of 
parameters of integrable cases, but also to find special asymptotic solutions even in non- 
integrable cases |^T|. The Painleve test finds wide use in theoretical physics as procedure 



for analysis of complex systems, for example, the self-dual Yang-Mills equation [£2| or the 
Bianchi IX cosmological model [53-56]. 

If one substitutes in a system of ODE's, for example, in system (2): 

X = ^a{t - to)" + E ^j+a{t - to)'+" and y = dy^{t - tof + J2 ^j+pi^ " toY^^, 

where N^ax, P, cx^ and cy^ are some known constants, then the system of ODE's 
transforms into a set of linear algebraic systems in coefficients cxk and c^^. In the general 
case one can obtain the exact solutions (in the form of formal Laurent series) only if 
one solves infinity number of systems: Nmax = oo. On the other hand, if one solves a 
finite number of systems then one obtains asymptotic solutions. With the help of some 
computer algebra system, for example, the system REDUCE [57-59], these systems can 
be solved step by step and asymptotic solutions can be automatically found with any 
accuracy. But previously one has to determine values of constants a, /5, cx^ and cy^ and 
to analyze systems with zero determinants. Such systems correspond to new arbitrary 
constants or have no solutions. Powers at which new arbitrary constants enter are called 
resonances. The Painleve test gives all information about possible dominant behaviors 
and resonances. Moreover, the results of the Painleve analysis prompt cases, in which it 
is useful to include into expansion terms with fractional powers of t — to- 

^They also were the first |3^ to point out the connection between the nonhnear partial differential 
equations (PDE's), which are soluble by the inverse scattering transform method, and ODE's with the 
Painleve property. They have proven that if a PDE is solvable by the inverse scattering transform and 
a system of ODE's is obtained from this PDE by the exact similarity reduction then the solutions (of 
this system of ODE's) associated with Gel'fand-Levitan-Marchenko equation will possess the Painleve 
property. Furthermore, they conjecture that, when all the ODE's obtained by exact similarity transforms 
from a given PDE have the Painleve property, perhaps after a change of variables, then PDE will be "in- 
tegrable" . Subsequently the Painleve property for PDE was defined and the corresponding Painleve test 
(the WTC procedure) was constructed by J.Weiss, M. Tabor and G.Carnevale |^ (see also [41-48]). 
For many integrable PDE's, for example, the Korteweg-de-Vries equation |42], the Backlund transforma- 
tions and the Lax representations result from the WTC procedure |4^, 44 1. Also, special solutions for 



certain nonintegrable PDE's were constructed using this algorithm |47|, tlsT The Painleve test for both 
ODE's and PDE's which is based on perturbation theory is presented in ||49|, IsT 
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The test consists of three levels, if the system passes some level then it means that 
the system may be a system of P-type and we have to check the system on the following 
level. 



3.2 The first level: find the dominant behavior 



We assume that the dominant behavior of solutions in a sufficiently small neighborhood 
of the singularity is algebraic. To find the dominant behavior we look for solutions in the 
form 

x = ai{t-to)"' and U = a'2{t - t^Y ■, (6) 

where to is arbitrary. Substitution (6) into (2) shows that for certain values of a and (3, 
two or more terms in the equations of (2) may balance (these terms have the same powers), 
and the rest can be ignored as t — > to- For each choice of a and (3 the terms which can 
balance are called the leading terms. For the Henon-Heiles system in the general form 
there exist two possible dominant behaviors |2|: 



Case 1: 


Case 2: {(3 < 3?e(a) < 0) 


a = -2, /3 = -2 
ai = ±3^2 + C, aa = -3 


a=^^v/i-^^/^ P= 2 
fli = arbitrary, ^2 = — ^ 



It is possible that an original system is not of P-type, but, after some change of 
variables, the obtained system is of P-typeQ. At C = — y in the Case 2 we obtain 
a = — I and to continue the Painleve test we have to introduce new variable z = and 
to consider the following system instead of system (2): 




1 ^2 



zi - 2Xz' - 4z'y, 



-y- z- 



(7) 



3.3 The second level: find the resonances 

The second level of the Painleve analysis is finding the resonances. For each obtained 
pair of values of a = 2a and (3 we construct the simplified system that retains only the 
leadings terms of equations of the original system (7). 
For a = —4 and /3 = —2 the simplified system is 



ZttZ 

ytt 



Substituting 



z 



54, 



to)-^ + bi{t ~ t, 



0) 



-4+r 



^Z^ 
2^t 



4z^y, 



and 



y 



-3{t - to)-"^ + b2{t - t 



0) 



-2+r 



^In this case the original system is said to have the weak Painleve property. For example, in the 
integrable case (iii) system (2) has the weak Painleve property. The interesting example of a system with 
the weak Painleve property is presented in BG] . 
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in system (8.1), we obtain that to leading order in (t — to) {(^"^o)'^ ^ ^oi the first equation 
and {t — toY~'^ for the second equation} this system reduces to a system of two hnear 



algebraic equations: 



Q{r)b = 0, 



(9) 



where Q{r) is 2x2 matrix, which elements depend on r, and b = (61,62)- Determinant 
det{Q{r)) is a polynomial of order 4. Equation (9) has nonzero solution only if 



det{Q{r)) = 0. 



(10) 



The roots of equation (10): 



ri 



r2 = 6, 



^^3 



5 + 753:8 



753:8 



determine resonances (one root is always (—1), it represents the arbitrariness of to). Some 
roots of (10) are not integer. This result means, that the Henon-Heiles system with 
C = — y is a nonintegrable system. There is no algorithm to find the general solution 
for a nonintegrable system. 

To find special asymptotic solutions let us consider the dominant behavior in the Case 
2 {a = —3 and (3 = —2). The simplified system is0 



Substituting (ai = a^) 

z = ai{t - toy^ + bi{t 



tr 



ZttZ 

ytt 



(8.21 



and y 



15, 



to)"' + h2{t - to) 



-2+r 



into (8.2) we repeat calculations and obtain that resonances and corresponding arbitrary 
parameters can arise in terms proportional to {t — to)^~') where r = —1, 0, 4, 6. Root 
r = — 1 corresponds to arbitrary parameter to, root r = corresponds to arbitrary 
parameter di, other roots correspond to new arbitrary parameters, i.e. new constants of 
integration. For arbitrary C values of resonances r are ^ : 



Case 1: 


Case 2: 


a = -4, 
/3 = -2, 

^ _ 1 fi 5 , Vl-24(1+C) 


a = ±.y/l - 48/C, 
/3 = -2, 

r = -1,0, 6, Tv/l -48/C. 



^The simplified systems (8.1) and (8.2) are different. 
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3.4 The third level: find the constants of integration 

The third (and the last) level of the Painleve test is a substitution into the original (not 
simplified) system (7) the following series: 

z = ai(i - to)"' + E c5j-3(t - io)'"' and y = - —(t - to)"' + E %-2(^ - ^o)^"', 
i=i * i=i 

where r^a^ = 6, cy^ and ci^ are unknown constants. 

After this substitution system (2) is transformed to sequence of systems of linear 
algebraic equations. Solving these systems we find cy^ and czj. Determinants of systems, 
which correspond to resonances, have to be zero. 

For example, to determine and czi we have obtained the following system: 

' 557056S? + 5f(l5552000A - 4860000) + 

< a? (8640000005^2 + 108000000A2 - 67500000A + 10546875) = 0, (11) 
818176at + a?(l5660000A - 4893750) - 8100000005^2 - 6328125 = 0. 

As one can see this system does not include c5i, so czi is an arbitrary parameter (a 
constant of integration). For any A this system can be solved as a system in 5^2 and Oi. 
We obtain new constant of integration c5i, but we must fix Si, so number of constants of 
integration is equal to 2. It is easy to verify that cy^ is an arbitrary parameter, because the 
corresponding system is equivalent to one linear equation. So, we obtain an asymptotic 
solution which depends on three parameters, namely to, czi and qj^. 



3.5 New asymptotic solutions 

Now it is easy to obtain asymptotic solutions with arbitrary accuracy. For given A one 
has to choose oi as one of the roots of system (11). After this the coefficients cZj and cy^ 
as functions of cz^ and cy^ can be found automatically due to computer algebra system 
REDUCE. For some values of A asymptotic solutions have been found as the following 
series (without loss the generality we can put to = 0): 

50 "l^g 50 

z = ait-^ + E czj_3t^-^ and y = - —t'^ + E %_2^^"^- 

For example, if A = ^ then system (11) has the following solutions (cii ^ 0): 

{ _ 25\/2 _ __ 1819 1 i _ 25iv^ _ __ 8700683 1 

16"' "" 663552/' " ' "" 1364926464] " 

If oi = ±^fg^ then we obtain: 

_ I 125^-2 I 25v^ .-l , 1625 i-p:Zf< 

^ 16 "^192'' 768 82944 "T t-^l^-T 

, ( 21845 _ V2~ \f2,( 437425^2 _ 25^2 ~ _ im_~ \ +3 , 
"""V 47775744 q ^■^l J ^ V 9172942848 48 ^^4: 3456 ^^1^'' "T • • • , 
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32 



205 I 115^2 y. 
2304 13824 



1819 
663552 



I ( 1673v^ , 1~ \ .3 I ^, ^4 , / 1044461^2 
^ V 11943936 ^ 6 ly '-^4'' ^ ^ 220150628352 



19 



and 



-2 _ 25v^ ^-l I 1625 I ~ 
768 82944 ""^1 

^ 48 



9216*^^1 2*^1/4) + ■ 



_ 25VI|-3 I 125^- 

16 ~'~ 192 768 

I ( 21845 I V2~ \f2_:( 437425v^ _ 25v^ ~ _ 191_;;r ^ +3 , 
"•"V 47775744 6 ^^1/ V 9172942848 48 3456^^1/ 



_ 15+-2 _ 5v^|-l _ 205 _ 
8 32 2304 

1673v^ , 1~ \ .3 : ~ .4 
11943936 ^ 6''^1J ^ ^^4'' 



115v^ 



t 



1819 

13824 " 663552 

/ 1044461^2 
I 220150628352 



19 
9216 



CZi 



23V2; 
384 



cy, + 



For real values of parameters and time these solutions are real. In the case ai = 
following table shows how cz^o and c^gg (coefficients of terms proportional to t^^ 
on the arbitrary parameters czi and cy^: 



(12.1) 



;i2.2) 



^the 

16 

) depend 



CZi 


cy^ 


5550 


cy50 


CZi 


cy^ 


5550 


ci/50 


-1 


-1 


4x10-12 


-1x10-13 








-1x10-^^ 


2x10-^5 


-1 


-0.6 


4x10-12 


5x10-1^ 





0.2 


-8x10-20 


-3x10-20 


-1 


-0.2 


-1x10-1'' 


-2x10-1^ 





0.4 


-2x10-1'^ 


-9xl0-i« 


-1 





-1x10-20 


3x10-22 





0.6 


-5x10-1*^ 


-5x10-16 


-1 


0.4 


6x10-1^ 


1x10-16 





0.8 


-5x10-15 


-2x10-15 


-1 


1 


6x10-12 


8x10-1^ 





1 


-3x10-1^ 


-1x10-1^ 


-0.6 


-1 


3x10-12 


-5x10-1^ 


0.4 





-4x10-25 


-1x10-26 


-0.6 


-0.6 


4x10-1^ 


-1x10-15 


0.4 


0.4 


-1x10-15 


2x10-1^ 


-0.6 





-5x10-23 


9x10-25 


0.4 


0.8 


-3x10-13 


2x10-15 


-0.6 


0.4 


3x10-1^ 


4x10"" 


0.8 





-1x10-21 


-3x10-23 


-0.6 


1 


3x10-12 


2x10-1^ 


0.8 


0.4 


-2x10-15 


2x10-16 


-0.2 


-1 


1x10-12 


-2x10-1^ 


0.8 


0.8 


-6x10-13 


1x10-1^ 





-1 


-3x10-1^ 


-1x10-1^ 


1 


1 


-4x10-12 


1x10-13 


20 


20 


-2.2 


0.051 


20 


40 


-603 


6.88 


40 


20 


-11.1 


0.01 


40 


40 


- 1128 


24.5 



One can see that coefficients tend to zero very rapidly when the absolute values of the 
parameters are less than unit. 



4 The connection between asymptotic solutions and 
exact solutions 



E.I. Timoshkova ||2J| found that solutions of the following equation 

2/2 = Ay' + ^2/5/2 + ^y2 ^ ^y3/2 ^ ^ 



(13) 



where A = V = 0, B, C, G and S are some depending on A constants, satisfy (3) at 
C = — After change of variables: y = g"^, we obtain the following equation: 

g^t = ^{Ag'' + gg^ + Bg^ + £g + C). (14) 

Equation (14) is of P-type. The general solution of equation (14) has only one arbitrary 
parameter to and can be expressed in elliptic functions |6^ . 

For all values of A the trajectories of motion are given by the following equation: 

^ ^^Gy'/' + 2{B + l)y + + c] . (15) 



= z = — 



2 V2 " ' 2 



Let us compare the asymptotic solutions (12.1) and (12.2) with these exact solutions. 
It is easy to verify, for example, with the help of the Painleve test of equation (14), that at 
points of singularities these solutions and our asymptotic solutions have the same asymp- 
totic behavior. It means that for some values of parameters czi and cy^ our asymptotic 
series give exact solutions and, hence, converge. 

Let us consider in detail the case of A = |. In this case the following solution has been 
found [0]: 

where y is a solution of the following equation: 

yl + §2/^ + iy^ ± -^y"^ = 0. (16) 
^* 15^ 9^ ^135 

Depending on a choice of a sign before the last term, we obtain either (in case of 
sign +): 

25(1- sin 

3(l-3sin(^))' 9(l-3sin(^ 
or (in case of sign — ): 

5 25fl + sinf^^' 
y = ^ and = ^ ^ ' ' \ . (17.2) 

3(l + 3sin(^)) 9(l + 3sin(^^)) 

It is easy to verify that the series (12.1) with 

858455 , 3205v/2 

cyA = and czi = , 

^ 12039487488 3981312 

are the Laurent series of (17.1) and the series (12.2) with 

858455 , 3205v^ 

cyA = and czi = , 

^ 12039487488 3981312' 

are the Laurent series of (17.2). So, our asymptotic series converge at these values of 
parameters. It is possible that our asymptotic series converge also for other values of 
parameters. I plan to analyze this question in future. 



and = ^ \ ; (17.1) 

3 
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CONCLUSIONS 



Using the Painleve analysis one can not only find integrable cases of dynamical systems, 
but also construct special asymptotic solutions even in nonintegrable cases. 

We have found the special solutions of the Henon-Heiles system with C = — ^ as 
formal Laurent series, depending on three parameters. For some values of two parameters 
the obtained solutions coincide with the known exact solutions. 

The author is grateful to R. I. Bogdanov and V. F. Edneral for valuable discussions 



and E. I. Timoshkova for comprehensive commentary of [24|. This work has been 



supported by the Russian Foundation for Basic Research under grants A/° 00-15-96560 
and M° 00-15-96577. 
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